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I.   BACKGROUND 


i  he  overall  research  o-f  which  thie  thesis  forms  a 
has  a  three— -fold   purpose. 


i  r  sci  v  ,   it 


tne  student:  about  me 


s amenta. 


intended   to 


govern  the  aerodynamics  of  wings   in 
jndly,   it   teaches   the   student 


bee: 


ncomorsssi b  ± e 


ow   to   convei 


fundamental  equations  into  a  form   suitable   for   numerical 


calculation  on  me  digital  computer,  ana  to  ns.rry  ou" 
complex  programming  required  for  this  purpose.  Th: 
' rst  two  obiectives  have  been   accomplished 


tne 


c  omo  u  t  er  or  os 


' o v i d e s  a  u setui   d e o a g o g i c  a  i 


desion  tool  for  illustrating  the  effects 


-ign  paramenters  on  the  final  aerodynamic  perform; 


a  wi  nq 


ucar 


.  owar 


J.  L.  Parks  made  a  creditable   start 
first  objective  in  his  thesis  CRef  .  U,  However ,   in 

limited  time  available,  he  was  unable  to   make   sigmf 
progress  toward  the  second  or  third  objectives. 

It  was  evident  early   in   this   present   analysis 
Parks'   computational   technique   had   some    flaws, 
expectation   was   that   these   flaws    would    be    re 
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ne 


a  d  l  i  y 


discovered  and  soon  corrected  and  that  objectives   two   and 
three  would  be  rather  ouickiv  and  easilv   comoleted. 


at 


nas  not  oeen  tne  case : 


H5  "cne  investiqation  actual  iy  uevsiopeG.  a  whole  series 
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£Ai  .— i  "— . 


encount  ef 


^nn       ee~n        waj 


evsi 


■  d  ii  n  u       inan 
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obstacle  was  really  attributable  to  some  conceptual  error. 
Once  the  error  was  "found,  the  obstacle  disappeared  and  the 
investigators  gained  a  deeper  insight.  Some  of  the  major 
obstacles  that  were  encountered  and  eventual Iv  overcome  are 


Confusion  anouTi  une  proper  type  of   computational   grid 
to  use,  whether  staggered  or    unstaggered ,  and  the  grid  size 


necessa 


t—  j_ 


-^  .->..'    +.-,-- 


or    accuracy  was  a   majo 


-ir-.~-     s 


actor,     «n   fci 


tor   oooc 


unstaggered  grid  over-    tne  semi —span  is  now  usee, 
reason  =. 

Another  point  of  con-fusion  was  how  to  deal  numerical  I  v 
with  the  singularity  that  occurs  in  the  governing  integral 
equation*  One  option  is  to  avoid  the  problem  by  staggering 
the  grid  o-f  field  points  with.  respect  to  the  end  of 
control  points.  The  other  option,  which  was  chosen  for 
this  investigation,  resolves  the  indeterminacy  by  rigorous 
analysis  and  as  a  consequence  allows  the  employment  of  a 
simple  unst aqgered  mesh. 

A  third  point  of  confusion  concerned  the  evaluation  of 
the   partial   derivatives   of   the   circulation    function. 


Pe   usee 


Should  analytical  or  finite  difference  method' 

represent  these  derivatives^"   It  was  found  that   analytical 


?ntiation  is  i  ncorrect  and   that   finite   aifferences 


iiiiii  anoiher  contusinq  ocint  cancarned  the  validity  o-r 


while  this  procedure  is  widely  advocated  in  the  technical 
literature,  it  was  -found  that  a  much  simpler  and  clearer 
formulation  can  be  obtained  otherwise* 

Boundary  conditions  also  caused  some  confusion. 
Special  conditions  apply  at  the  leading  B.nd  trailing  edges 
and  at  the  wing  tip  and  midspan.  These  quite  complex 
boundary  conditions  have  now  been  fully  and  rigorously 
analyzed  and  incorporated  into  the  formulation  of  ihe 
or ob 1  em ■ 


evident  from  the  foreqoinq  discussion 


research   nas   amounted 


Tin 


GGLicat  i  or 


aer odynamics  and  in  numerical  methods.   In   these   respects 
it  has  been  a  richly  rewarding  experience. 

However.  in  view  of  the  foregoing  obstacles  and 
problems,  the  time  schedule  has  of  course  been  greatly 
delayed  from  that  which  was  initially  anticipated.  Thus, 
there  a^e  no  final  numerical  results  at  this  particular- 
stage  of  the  investigation.  There  is  also  no  more  time 
available  to  this  investigator.  Hence,  this  final  aspect 
will  have  to  be  completed  by  some  subsequent  investigator. 
Nevertheless,  what  the  present  effort  has  produced,  to  aid 
any  subsequent  worker,  is  a  very  soohi st i cated  and   refined 
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strength  over  the  wing,  trailing  to  infinity  aft  of  the 
wing  is  utilized  to  determine  the  value  of  the   circulation 

function  at  a  finite  number  of  control  points  on  the  wing. 
The  vortex  sheet  strength  is  restricted  by  the  Kutta 
condition  at  the  trailing  edge  of  the  wing  and  by  the 
requirement  of  no  flow  throuqh  the  wing   at   the   specified 


conn 


r est r  i  ccs! 


steady,  inviscid,  incompressible  now  about  tmn  wings  witii 
straight  or  swept  leading  and  trailing  edges. 

Wing  aerodynamics   involves   two   fundamental   aspects. 
First,  a  desired  basic  lift  Distribution  mav   be   specified 


over  the  wing  and  tne  corresponding   wing   camber   must   be 
determined.   Second,  a  flat  plate   wing   may   be   specified 

with  the   resultant   additional   lift   distribution   to   be 


c a i cui ated 


Jnce   the   additional   lift   distribution 


known,  other  pertinent  aerodynamic  parameters  may  di 
determined;  notably,  the  slope  of  the  wing  lift  czu.r^-^  B.m 
»■*■  the  aerodynamic  center 


tne  location  o- 


hi so ,  wnen  tne 
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chordwise    ores5u.r?    d i  s t r  i  b u 1 1  o n s    q ver    ths 


wmch    srs    necessary    in    desioninq    the    wma 


.iciur 


Althouah   Parks'   calculations   did    not 


actory  fssuit ,  he  pointed  out  areas   of   concern 


oossipie   errors   ana   maicatef 


need 


Ulhere    Harks       used       a       staggered       mesh       between       o: 
points       and       -field       points       on       the       wine       to  avoid 


"■aerer mi  nac  v    encountered       whenever 


H'-' 


rieia  point:  coincide,  tne   present   analysis   reso. 
indeterminacy  and   works   with   control   points   a? 


points  superimposed  upon  each  other.  Also.  the  series 
solution  employed  by  Parks  to  represent  the  circulation 
distribution  has  been  discarded.  Instead,  this: 
distribution  is  now  represented  more  simply  and  di recti v  bv 
the   values   of   the   circulation   at   the    fielo    points 


"c  n  e  m  s  s  i  v  e  s 
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cqr5I5ls  of  a  wing  plan-Form  symmetric   aoout   miaspan   with 

straight   leading   and    trailing    edges    containing    no 
discontinuities,  except  at  midspan  where  allowance  is   made 


or  sweepback 


ihe   winq   nas   no   control   surfaces. 


H 


typical  wing  is  pictured   in   rig.  1.   with   tne   geometr: 
coordinates  labeled.   The  winq  is  completely   specified   [ 


eaci no 


e  o  o  s    s  w  e  e  d 


nates   on   tne   winq   sr^       normalizes 


semi— scan  lenqtn  ot  one  uni 


Sweepback  at  the 


Lamda(A).  Eta(T^)  is  the  normalized  spanwise  coordinate. 
Parameter  Tau  (f)  defines  wing  taper.  Xi  (1  )  is  the 
normalized  chordwise  distance  -from  the  leading  edge. 
Sigma  (CD  is  a  sign  parameter,  +1  -for  the  right  wing  and  —  i 
-for  the  left  wing.  The  mean  chord  (c)  is  a  function  of 
aspect  ratio(AR) ,  that  is,  c  =  2/AR.  This  analysis 
involves  two  coordinate  transformations, 
transformations  proceed  from  the  (x,y)— plane  to  tni 
non-dimensional  (|,\'-plane,  then  to  the  i<p,6) 
coordinate  Diane. 


ne 


incui  a 
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Figure    1 .      Wing  Coordinates 
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'_!__'!    '    U  < 
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klp  =    h\? 

c    =     C2/AR)     C  (  1     +T)     -    2Ttrn.] 
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(2/AR)2T<rdTy- 
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cdl    +    Idc 


=    Ao"dr\_  +    cd\    +    |  L  -  (2/  AK)  2T<rdy\_  J 

yU  =    C   A  -  C2/AR)2/rl  j 

dx     =    cdi    +y^CTaT\_ 

T   =     (i  —  taper  ratio) /(l  +  taper  ratio) 

This   transformation   simplifies   the   algebra 
problem  significantly.   The  chordwise  dimension  now 
•from  0  at  the  leading  edge  to   +1   at   the   trailing 
while  the  spanwise  dimension  varies  from   —1   at   the 
wing  tip  to  + 1  at  the  right  wing  tip. 
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*m  tns 
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;sing  tns  relationships 


c    =    0. 5    *     ( 1    —    cos  0  )  ;        d| 
,  \   =    cos  8  ;        d\=    -sin   6     dB 


si  n  <2>    d0 


1  6  J 


rom    eauation     'o=,  li>;  „     a 


soqe    ro     i     at    trie    irauina    eoas 


s     \    varies    f 

0    wi  1 1 


■3.  1 


m    zero    at    the       i  e a ti i  r~' q 


v  a  r  v 


-from    zero    to  TT\ 


espsct  i  vel  y .        Froin    equation     (3.  17)  ,     as    V\  varies    fr 

0  w  i 


se  i  e- 


w  i  n  q  1 1 


at  tne  nqnt  winq  t 


t'cm  TT  to  zero,  respectively.  These  limits  Bt~^  important 
in  subsequently  developed  integral  equations.,  An  added 
bonus  of  this  trans-Formation  is  that  it  results  in  mesn 
points  near  the  leading  and  trailing  edges  ana  wing  tips  to 
be  closely   spaced,   allowing   -finer   resolution   in   tnese 

-ig.  2  shows  the 
coordinates   (£*YT> 
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Figure  2.   Linear  Coordinates  ( ?  ,^l)  versus 
Angular  Coordinates  (0,6) 
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IV.    bOvLKNINb  tLiUAiiGNb 


The  vertical  velocity  induced  at  any  particular  control 
point  consists  o-f  that  contribution  made  by  the  vorticity 
distribution  over  the  wing  plan-form  plus  the  contribution 
made  by  the  trailing  vortices  behind  the  wing.  Parks 
CRef  ,  1]  provides  a  quite  detailed  derivation  o-f  tne  two 
general  equations  for  induced  velocity 


i-"_p  l  i 


?  o  u  a  tions  m  a  v 


also  oe  round  in  many  Diner  aeronautical  references  such  as 


Bertin  and  Smith  CRef.  k'j  and   Kuethe  sr 


now   Li-\ej 


"he  fallowing  form  of  these  equations  was  utilized  in 


anal ysi s 
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jp 


^JJ  r^Mfr^r^^ 


where 


4rrJ    (^-fp) 
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3V.a1 


r~    =     i.x     —    x     )  ■*■    +     (y    —    yryijC" 
and    rj.'    is    obtained    by    replacing       x       with       >-.  r 


PHM^Tl r 


(4,  3)  . 

Equation  (4.1)  is  basically  the  same  as  Park's  equation 
(3.6)  where  w  '  is  the  velocity  induced  at  a  control  ooint 
due  to  ~th<s  vorticity  distribution  over  the  wing  surface. 
Equation  (4,2)  is  the  same  as  Park's  equation  (3.8)  where 
w"  is  the  velocity  induced  at  a  control  point  due  to  the 
trail inq  vorticitv  behind  the   winq,    Subscript    t    here 
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A,   CHANGh.  OF  I"7  TO  WING  COORDINATES 

Circulation   has   the   dimensions   o-f   velocity    times 
distance.   Consider  two  separate  circulation  -functions  =  one 


ion — cJi (oensx onai  with  respect  to  semi- sna.n  and   tree   stream 


ve i oc lty.  ana  one   n  on -a i mens i on  a i   wit; 
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i^s  r(x,y)  =  v^c  r(f,^)5   E/s  =  2/AR 
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Equation  (4, 1)  may  now  be  written  in  the  following  form 
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B.   CHANGE  OF  U   TO  ANGULAR  COORDINATES 

Utilizing  equations  (3.16)  and  (3.17)  together  with  a 
non— di  rnensi  onal  circulation  function  in  the  \(p=6)—  plane  ti 
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Upon  reversing  the  limits  of  integration  and  reversing  the 
algebraic  sign  of  the  integrand,  equation  (4.23)  can  be 
wr i tten  as 
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-^-nP)(g))d^da 
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tquatiDn     (4.^:4/     is    the    governing    equatioi 
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velocity  at  eacn  aanzroi  point  nu&  ig  me  summation  of  the 
effects  of  every  field  point  on  the  wing  plan-form.  Note- 
that  all  factors  inside  the  integral  ar^  a  function  of  wing 
olanform  geometry  exceot  for  the  two  cartial  derivatives  of 


:ne  circulation  i- unction,   hi  so  note  "c; 


becomes   indeterminate   when   x  =  ;■-'    and 
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coincide  and  is  discussed  latere 
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the  trailing  edge  and  simply  require  the  subscript  "t'  to 
be  added  to  the  appropriate  variable,  However,  since  £  is 
constant  at  +1  along  the  trailing  edge  a  few  of  the 
equations  change  slightly  and  sn=r  rewritten  here  for 
cl  an  ty 
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(3.12)  becomes  dxt  =  l  A  -  (2/ AR)  2T3(Tdr^  .(4=27) 

( 3 .  14)  becomes  d:u  =  AO"dn  ( 4 .  28 ) 

The  circulation  along  the  trailing  edge  of  the  wing   is 

a  function  of  the  spanwise   coordinate   Y^   only. 
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d6 


v^t-,  o-t>  > 


Equation     (4.35)     now    becomes 


a\   _i_  as 


wp     =  4ir\Aft/4  (n-V>    ^9 


d£ 


(4. 38) 


Aqain,  reversing  the  limits  of  integration  and  the 


e.i  qeorai  c  si 


gn  of  the  integrand,  the  integral  equatis; 
velocity  due  to  the  trailing  vorticity  be? 


.IT 


w. 


' /-*) 


£_  m)  de 


-4rr  Wv/  K' V)  IdS 


4.  v 


D.   TOTAL  VORTICITY  EFFECT 

The  total  induced   velocity   at   = 
noint  is  the   summation   of   the   effect   produce 
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oy 


vorticity  over  the  wing  and   the   effect   produced   by 
trailing  vorticity.    This   means   that   for   each   cont 
point ;  the  right  hand   side   of   eguation   (4,24)   must 
evaluated  at  every  field  point  ovsr    the  entire  wing  surf 
and  summed-   This  sum  will  give  the  effect  produced  oy 
vorticity  over  the  wing.   Also,  for  each  control  poini. 
slight  hand  side  of  eguation  (4,39)   must   De   evaluated 
each  station  along  the  entire  trailing  edge  of  the  wing 
summed.   This  sum  will  aive   the   effect   produced   bv 
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evaluated,   This  may   be   expressed   as   a   combination 

equations  (4,24)  and     (4,39)  as  follows,   where   w,_   is 
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pianform;  however ,  symmetry   about   micspan   again 


s  t=i  >-—  ;Ti : 


utilization  o+  numerical  values  over   only  an-^       semi -spa 
The  nqht  semi— soan  is  chosen  for  this  analysis, 
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mesh  coincide,  Equation  (4.39)  also  becomes  indeterminate 
when  the  spanwise  coordinate  of  the  trailing  vorticitv  and 
the  concerned  control  point  coincioe.  Note  that  G  in 
equation  (4,34)  approaches  zero  as  the  spanwise  coordinate 
of  the  control  point  and  the  trailing  vorticitv  converge. - 
One  method  to  avoid  these  indeterminate  points  is  to 
staogsr  the  different  meshes  so  that  they  n^v^t~  coincide. 
Another  method  is  to  determine  an  analytical  limit  or 
principal  value  of  the  integral  as  the  indeterminate  point 
is  approached.  The  second  method  is  chosen  for  this 
analysis  and  is  presented  in  the  following  paragraphs. 
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where  A  and  B  are  known  -functions  of  A?  0  n  ?  0ri  ^  ^^  *.  fD 
and  AR.  The  denominator  o-f  the  integrand  o-f  equ< 
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Assume,  by  analogy  with  the  above  development,  that  the 
numerator  o-f  equation  (4=24)  may  also  be  expressed  as  a 
-function  o-f  u  and  v,  namelv 


F   =  Du  +  Ev 


(5. 20) 


Now,  switching  to  polar  coordinates  with   /»  '  the   polar 
radius  and  '<*   the  angle,  let  u  =  P  cos  04  and  v  =  o    sin  oL    = 
The   e-f-fect  o-f       the   singular   -field    element    must    be 
determined =    An   integration   over   the    entire    element 
produces  the  reauired  result.    The   intearation   covers   a 
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where    H  (.<*■)     is    also    an    abbreviated    form. 
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polar   symmetry.    This   proves   that   the  effect   of   the 

singular  field  element  on  the  control  point  i  <. 
be  neglected. 
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As  t\   approaches  ri_.  in   equation   (4=39) 


IDpi 


'  f=  f  r~ir 


the  inteqrand  is  indeterminate.   Note  that 


foes   nor   approach   zero,   since 


:he 


x-cooroinate  of  me  center  ot  tne  concerned   control   ooin- 


and  H+.    is  the  x— coordinate  at   thi 


aiiino   edae   ot   the 


wi 


i  iy 


n   the   immediate   vicinity   ot   :ns   sinau. 


'iTidi  ns 


( V^  —  T^n )  becomes   very   very   small,   but 
therefore,  it  can  be  represented  as  a  smal i  value   timf 
constant.   Let  this  constant  be   (x±  —  >'n);   therefore 
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Thus,  in  the  limit,  as  ~Tl  approaches  Tl  ,  e  approaches 
zero,  and  the  value  of  equation  (5.32)  becomes  zero.  This 
proves  that  the  trailing  vorticity  contribution  directly 
downstream  of  a  control  ooint  makes  no  contribution  to   the 
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This  analysis  involves  mapping  points  on  the  wing 
the  physical  (x,y)  — plane  to  corresponding  points  in  the 
(0,0)— plane  and  vice  versa,  Fig,  4  depicts  a  uniform 
rectangular  grid  o-f  points  in  the  ((p,&)  -plane,   constructed 
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Figure  4-.   Computational  Grid  in  (0,6)  Plane 
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9,  =  Ae^J-0.5);   A8  =  2cT  =TT/I6: 
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The  above  equations  will  produce  64  grid  points  on  the 
right  semi -span  with  each  point  centered  in  its  respective 
rectangular  mesh  block.  The  ' j  subscript  varies  spanwise 
from  right  wing  tip  to  midspan  while  the  i '  subscript 
varies  chordwise  from  leading  edge  to  trailing  edge,  The 
unknown  values  of  the  circulation  function  at  each  of  -chess 
g r  i  d  p o i  n t s  must  b e  deter m i  n ed , 

For  every  field  element  on  the  right  wing  there  exists 
a  corresponding  mirror  image  field  element  on  the  left 
wing,  By  utilizing  the  symmetry  between  each  field  element 
and  its  mirror  image,  the  range  of  the  ' j '  subscript  can  be 
restricted  to  8  vice  the  16  it  would  otherwise  be.  Fig,  5 
depicts  the  plan-form  distribution  of  the  circulation 
function   with   appropriate   subscripts   for    matrix    and 
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vice  versa  the  following  relationships  a.r^    used,   with 
denoting  the  respective  vector  subscript. 
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Figure  5«   Gamma  Matrix  and  Vector  Grids 
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Pv(64-) 

CD 

UJ 
en 


-o 

UJ 


u 
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k/bS    bv    intecer    division 


where  IM< 


)   =   P  (fc) 


(6,6. 


(6.7) 


VARIABLE  NOTATION 


HQI.3I.  1  DH     1  D  Tl  r  DO  LI  C  ?  G 


+-hf=> 


Daraqraphs  is  used  extensively  in  the  subsequent 


The   following 

d^r i vat i ves  is  al so  used 


notation    tor    various    functions 


Note  the   different   exoressions 


correspoi 


:ne  same  numerical  value, 


The  first  derivatives  of  the  circuli 


men  en 


zv 


=   PAu,i>  =  V,a 


Z&l       0 


0 


($)  -  *«:..»  -  r.¥« 


cr        (xi     3  •_       '_j  f=:  i     a 


vat i ves  at  the  I eadi  ng  and  trai 1 i ng  edge 


x  6  =  1  U  > 


(6,11 


!  he  second  derivatives  at  the  I  eadi  no  and  trai  lino  edae 


J'P 


=  ft 


(  1  ) 


^     '0^_'J 


(  i) 


/  i.   1  71 
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The    lift    distribution 


lesiqnated    as    -follows 


AC,  ? 


I,    T 


r\.iki 


(6= 14; 


I  ne    I 


■  a    a  n  o    tra.iiina    e  a  a  e 


A£Vlj     =    PL(J 


An    )  . 


C'         {    -i  \ 


(6, 15) 


6.  16) 


C.   LIFT  DISTRIBUTION 

Consider  Bernoul  I  i  's  eqijation  tor  the  pressure  above 
the  wing  and  below  the  winq 


+     (C/2) 


I 


<oo 


4/2- 


=     Pi 


^/2>CVfi 


flti  /23" 


?6. i  / ; 


Rearr anqi 


-L    I  5  \_3         QHlj 


(Pl-Po)     4p 


J-    v/2 


T> 


=  A 


V- 


YL 


tor 


A. 


_2_  j_   ltfy 
MJ~   ah  c    Ulj 


i.  6 ,  i 


and  , 


an  =      2     ^n' 


(6.20) 


Now,  combining  terms  -from  above 


AC, 


2^)' 


2.   \ten\_  8 


lAR/C   \Sin.0/U#/~    AR    C   151110 


in 


Or,  in  the  natation  lntrDGucea  in  section  b 


P(_i,i)  =   -=-  \<   *^ : >  (6.22) 


D.   LEADING  AND  TRAILING  EDGE  REQUIREMENTS 

The  boundary  conditions  pertaining  to  equation  (6.22) 
■for  the  basic  or  ideal  lift,  and  the  additional  lift, 
determine  the  nature  o-f  the  circulation  function  and  its 
derivatives.  The  kutta  condition  -fixes  the  pressure 
coefficient  at  zero  at  the  trailing  edge-  For  the 
additional  lift,  the  pressure  coefficient  is  infinite  ai 
the  leading  edge.  However,  for  the  basic  lift,  the 
pressure  coefficient  at  the  leading  edge  is  finite*  These 
spec i f i cat i ons  together  with  equation  (6.21)  or  (6.22) 
determine  the  nature  of  the  circulation  function. 

First,  consider  the  Kutta  condition  which  requires  that 
the  pressure  coef+icient  be  zero  at  the  trailing  edge.   The 


ion  (6.22)  is  ai so  zero  at  ine  trai 


in+inite, 


,  since  <p  =  TT  and  consequently,  =^\r,  <p  —  0.  This  means 
I  f£  must  not  be  finite  or  the  quotient  becomes 
Therefore,  IJ^  at  the  trailing  edge  must  be  zero, 
resulting  in  an  indeterminacy.  Using  L'Hopital 's  rule  co 
determine  the  limit,  which  must  be  zero,  the  denominatior 
becomes  cos  Tf  =  —  1  and  the  numerator  must  be  zero.  This 
fixes  the  value  of  the  second  derivative  of  the  circulation 
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unction  with  respect  to 


^  a- 


Z  f  f   D     •=•.  i  S  O  i 


;.  u  m  m  — ;  r 


-rorm 


Ac, 


i  !)     =    u 


^2 


J  >    =   u 


GB 


the       condition  tha.t  the 


rnpf  t 


A  *_  A  CTI  !  *_       A 


ai   the   iea 


U  i  !  Si 


additional  I  i  f  t ,   At  the  leading  edge,  <p   =  0  and  sin  0  = 
there-fare,  the  denominator   of   equation   (6.22)   is   ag; 

Thus,  l/ij  is  finite   at   the   leading   eoge   f or       1 


^  -^r  '_■ 


additional  lift,   Summarized  in  equation  form 
AC  )  j   =  P,  ( j )  =  +  **> 


\a .  ^iz- 


in 


■fV 


.  ( j  >  >  o 


(.6-2/ 


Finally,  for  the  basic  lift  the  pressure  coefficient 


'.  o  t  n  a  t  ■*"  o  r 


the  trailing  edge  using  L'Hopita 


rule   results   k 


iecond  derivative  o- 


ie  circulation  function  witn 


to  <p   being  a  finite  value.   Summarized  in  equation  for 


Ac 


D  ' 


rLvj; 


iR) 


L/Cj 


(iflL-'V"'- 
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t=        HlDSPAh:    hNC    wING     flP    KtQJ  3CRLMEN  i  S 


hi  on; 


boundary 


Denoted  Dy  sudsi 


nr 


rRd)  =  o 


ah 


-  rR<i>  =o 


(iS,  .:.  I  ) 


At  mi  dspan ,  denoted  bv  suoscnpt  'm  ,  allowance  is   made 


or 


iisconnnuitv  in  si  one  o+  tne 


rem  ati  or 


unct  i  on 


Note  that  the  circulation  function  is  continuous  at 
midspan,  but  its  derivative  is  discontinuous  at  this 
I oc at i on , 


F=   BASIL  AND  AuDIilUNAL  Ul RlULA ! IuN 

the   concepts 


Hf 


unaerscanai no   o- 


UC.=M  l_ 


a  n  o 


aaax tional  circulation  is  necessa^v  in  tne  ana 


a   d  i  v  e  r 


his  n\  a  v   n  - 


prooiesn.  i  nese  concepts  ^rs    orietiy  revieweo  nere; 

Consider  the  circulation  function  I  (0,0 >  of 
wing  at  some  arbitrary  angle  of  attack  °<»  = 
expressed  in  terms  of  a  "basic"  circulation  -function 
\p\<f)..&)  and  a  "specific"  additional  circulation  function 
\/\<<P,Q)     35  follows 


n^>,<9)  =  rv0,e>  +  (oc  -  ^p)Vr,i4>,6> 


':   C?  =  -J-jL.  . 


where  («i  —  °^g)  is  termed  the  additional  angle  of  attack 
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It  is  usetui  to  desianate  tne   proauc 
■  s    the  "net"  additional  circulation  a.nti 
luantitv  Vr,i<t>-.d)        itself   as   the 

C'i  "f  £  8r9n 1 7.  at  i  no  equation  (6.  32)  w 


np=i  nn. 


5  n  e  c  1  T  "5 


;5DBCt     no    o<    gives 


f-cP<4,0> 

dot 


rfl«M: 


not  ice    tnat 


UipJ)  has  the   characti 


per  radian,  arid    that  it  is  numerical  Iv  equal  to   the 


wmcn       tne      nei 


A    I       1_  t_J 


iiation        would  reaczh 


extrapolated    to       an       additional       anqle 


raa 


J.   !    1  1 


ac 


: n e  aQditlonai  ana  i  e  o t   atiai 


is   restricted   to   small   values;    therefore,    the    net 
additional  circulation  remains  small  in  comparison  with  the 


specific  additional  circular. ion. 

All  other  aerodynamic  parameters  which  s.r^  proportional 
to  circulation  follow  a  linear  law  analogous  to  that  state? 
in  equations  (6=32)  and  (6=33),  For  example,  the  wing  lif- 
coefficient  Ci   conforms  to  the  expressions 


CL    =    CLB    +     <0*  _  °<B)CLf 


t  H  ,   . \±f  I 


an 


dC_L 

a* 


U 


f  fy 


nqain       note  tna- 


me 


secific  addition a 


ccstticient    L.  r,    is    icenticai       wi 


the       iitt       curve       sloce 


•  Us  /aeQ     and    i  =■    numerically    e 


to    me    vaiu.e    "c; 


is    n?i 


s.aaitiona.1     11 


,riri-.  ¥  -; 


jauld  reach  if  linearly  extrapolated  to 


additional  angle  o+  attack   of   one   raaian 
however  -   the 


C\  '_J  Cl  J.    !    I 


idditional        anqle       of        attack         (  o<  —  °^tu) 


resrri  cue*: 


e m a i  i  vsi  ub* 


I  t  l_J    C\  *__'  '-J  i. 


i  onai 


x  j.    '    \-       ai  c   ai  ^u 


:usseo  in  r;e- 


44 


II.         Bfcf>iERHL_    rKUtii-tr;    HiVHLrbib 


A.   PROBLEM  SPECIFICATION 

As  mentiorsed  prsviousiy,  the  calculation   sequence 
involve   determining   the   wing   slope   function    fro! 
specified  basic  lift  distribution,  or  determining  the 


i  dri  hut1!  nn  -f  r- n 


inuuion  trcm  a   speciriea   wi 


si  one 


1 nteres 


(_   i  i  cr  i 


■e   botn   the 


and   additional 


ai str  i nut i ons, 


The  values  on  the  right  wmq  scan 


ai  i 


calculations.   The  right  wing  span  is   divided 
chordwise   and   sight   spanwi se   stations,    r 
sixty— four  control  points  a-v-<=r    the  semi— span, 
the  condition  of  no  flow   through   the   wina 


l  ;-_si_i-_fcrLj 


'eiocity  musx:  balance  "Che  normal  componer 


~^e 


tree    stream    velocity* 


w. 


+    ^   sin   oc    =    IzAi 


betting   the   tree   stream   velocit 
approximating  sin  ot  with  <^,  for 


one 


i-i !  ■  i    !_ 


Wz/ox)     —  *<    =    w      =    w 
w  ™p         "p 


By    analogy    with    equation     (6.32! 
jnction    at    each    control     ooint    may 


se    overall     wine    si  one 


'eoreeenieG       as 
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L  Wi 


'_•__■!»!!_=>_•!  ?C1  )  ' 


!l  Tin 


p   +   (Ol  ~  0^R>  W. 


I-1     un 


Pti 


rs+iects  the  5I  gd?  of 


■_  1  ■  ■=         niiiy 


according  to 


r  si  ai 


^z/ix)     - 


eat  _ 


wh  lie 


w 


OH 


a  u  ■_!  1    '_j.  *__: !  i  -r\  i 


c.  u  >_  c\  t_  r:. 


Note       that       the       riant 


consi sts    at 


> "_'.  Ill         *_?    !  ^T  •_•  "_1  C%   '_   i    L_"  I    i 


39) 


solving  for  the   additional   lift,   the   only   unki 


a    -?£.  ■* 


24)  and  <4.39) 


:ne   Derivatives 


IJLCLlUji    1  •_*;  I   ! ■ 


i  hese  derivatives  5.r^ 


1  n 


usrfnb    o  1     c n e    1  n  1  c  1  a  1  1  v    u n Known    v a  j.  ubs 
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TOi 


K  n  o  w  n  v  a  1  u  e  s 


.Lti  ci  l  i.  u;  1 


_=  u  e\  •_  i 


ons  are  require 


?  r  a  s  o  j.  Li  c  1  o  n 


■_  1— ■  1 1  \_  1  «_?  j 


generates  sr-        equation   involving   terms   of   ~tn^ 
circulation    distribution.      This    produces    si: 
equations  in   sixty— four   unknowns   which   must   be 
simultaneously  to  determine  the  circulation   values, 
solution  involves  the  inversion   of   a   64   bv   64 
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•f^n     rp  i 


wnen    sdi vmo    Tins 


di  stri but i on  = 


ie       LiPKnawn 


se       wine       5ioci 


•h    Dv    a*f       man 


'  '_l  t_  A  *  ■  V  Sf  I 


the 


: *_|  :_<  i  F    E=JU 


f  ■  i_J  f""i  ulj ;  » i  i  j  T  i— :  l_         ;-.CQ  T  l  JA   i     i   I  i  i_  i_-       f~i  CT       i~*  T  i-_:  i   •  i   i !     f\*1    T  ;   •"  K  i 


In  developing  -finite  difference  equations,  a  polynomi 


throuqn  a  number  o+ 


n  t  s  in 


The  polynomial  used  in  each  case  must 


i  iruLii!  *=;'_; 


exist,  and  it  must  satisfy  the   boundar 


i  ons 


inite   difference   eouations   ut 


(_'.'_  j.  i  a  i.  =■ '_: 


th« 


!~  pnr  r 


difference  formulas  require  a  point  on  either  side   of 
point  in  question.   Hence,  the  central  difference  aqua' 
3tr^    valid   for   all   mesh   points   except   those   that 


lacent  co   che   leaomq   ana 


For 


i  e  a  a  i  n  o  e  u  g  e 


r  =  h,4>  +  ii/2)H-,4>2  + 


ore-"  j 


l~V  =  Al  +  A2^  +  DC^2:3 
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Equations  (7=6)  and  (7.7)  satisfy  ths  boundary  conditions 

-for  ths  additional  lift.    These  ars  that  P  vanish  and  l\ 

iere  <p    —    0-   Applying  equation 


d>  is 


t  p  £=      !?=.•=; 


/     ai 


qno    points    adjacent:    to    tne    i  ea 


sn   -£r-s.-   di 


for  Hi  anc  h- 


suhstitutinq  into  equati or 


tor  P  a(  in  terms  of  I  1  and  l  -3 


0  iP 


is  obtained  as  follows 


P^  =  a/A^)rP1  +  (l/3>r2l 


"herefore,  for  any  mesh  point  adjacent  to  the  Is 
:he  following  subscript  notation  is  used  for  )A 

fy<j,i)    =    ( i/A<£>  c  P<j  ,  i  >    +    <  i/3)r<j  ,2) : 


!,  /  , 


h  s  1  m  1 1  a 


-Poll  owi  no 


■/si  5  aopl  1  eo  7:0 


.  .■  t  cr   '_•  ^*.  =2  . 


r^(j,l)  =  <1/A^)  C-(A</>/2)2^L(j)  +  6p<j?i>3 


The  procedure  for  determi ni ng  the  difference  formula 
applicable  near  the  trailing  edge  requires  defining  a 
function  for  the  circulation  that  is  numerically  equal  to 
[  +.  at  the  trailing  edge  with  first  and   second   derivatives 

This   requires   beginning   with   a.   cubic 


equal  to   zero 


equation;  therefore,  1  e"i 


1=1.  +  (l/6)ATj(Tr  -  0 


)  ■-' 


11) 


Differentiating  once, 


48 


fy  =  -(i/2)H3^ir-^»2 


i srentiatina    aoai 


fytf    =    --3^-* 


iODi  VI  PG     BGLia" 


1  ) 


ira: i inq 


Bsoress: 


or  At,  and  substituting  into  (7.  12)  , 
■n  -for  I  ri    in  terms  of  I  (j,7)and  1  Cj.S 


i 


ait  t scenes 


i=  ±  ~  i 


i  QW5 


Fa<j,S)  =  (1/A0)  C  (-3/13)  i!j  ,7)  +  (3/13)  P<j 


R)  1 


Hi  SO 


rtu» 


=  (-1/26) li j ,7)  +  (27/26) l(j,S) 


ror  mesn  points  >'  through  /     in  the  cnordwi  ss  directio! 


LEnir3i  di  i  tG^encss  sr~s  usee  wn i  c n  resu. 


u  =>   J. 


n  me  toll  owing 


r^(j5i)     =     <1 /A<£)  C  (-1/2)  r<j  ?i-l)     +     Cl/2>P(j,i+l>  3     (7,16 


The    di "f -f erence    eouations    for    the    trai  1  ina    prioe    and    the 


; —  c^  r-v  -r-  y-  .3  ] 


i ,  ->  1  ■;  ,-i 


BPANWISE  DERIVATIVES  OF  CIRCULATION 

Aiona  the  riqht  wma   tin.   the   circula 


\_  i  «__'! 


a.r^    precisely  analogous  to  the   boundary   condit 
the  leading  edge  for  the  additional  lift* 
difference  sanation   develooed   above   -For 


ppce.     "cne       same 


5ese       Dounoar' 


w 1  nc 


Pft(l,i)     =    (i/A8)cP<i,i)     +    Cl/3)P(2,i)3 


ra^J?i>     =     (1/A0) C (-1/2)  P(j-l,i)     +     (l/2)rr( 


0 


1  , 


At    midsoan,     allowance    is    made?    -far       a       discontinuitv 


1DP         0  + 


L_  '_^  i    •=».   :_   i    !_■ !    » 


T u nction 


tTfEHTl  I  DOS? 


PV]  O L? ! 


Ha  a: 


poi  ynoflii  a j.     suitaoie 


+-  f-s         =  =>  *•  i  =  ■£  ^ 


at  i  stv       t ha 


(T-To)  =A,(e-e0>  +  <i/2)A^<e  -  ©o>2 


7 .  19) 


iQ    =    Aj    +    A  2  (  &  ~  O 


tsatis+yincj    ?o l.'. atiD n     v 


=    6    and 


50|  VI 


ri-{         aDG       H^ 


ittici  ents 


lv 


0) 


sDDivinq       eouaiio! 


inicspan    and    Equation 


20 ) 


a    L  •  I !  A    L_l  =•  !_'  C^.  I    1 


u"  C       o.  *C        1       —      ■_'      Q  X  V  3  zi> 


rm 


T(6,i)    -    (5/4>r<7,i 


)     +    (15/8)1    ( 


>r\a,i 


rem<i> 


C 1 /AG)  111(6,1 )     -    3T(7,i)     +    21(8,1)3 


Iq(8,D     =     (1/A©)  C  (1/2)  P(6,i  )-2r<7,i)  +  (3/2)  P(8,i)  1  (7.23: 
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t-iEl_D  ELEMENT  AND  IMAGE  ELEMENT 


Igebraic  signs  or  T^,   T, 


9 


and  T. 


Dart=,  one  nart   consist ina   of   tne 


iiriDuiiGr 


_\  ^.  e»  i    j 


ieiQ  Doint  on  tne  richt  semi— scan 


c-ntral  point,  and  another  Dart  consisting 


contribution   that   the   image   field   pom' 


;he   same   control   coin- 


concent  applies  to  ^  '! 


:he   trailing   vor ti city   Integra 
Ed  nations  (4=24 )  and   ( 4  =  3  9 )   sn 
modified  sliqhtlv  so  that  the  effect  of  a   field   point   oi 


aionp  ine  traxunq  e a d e 


tne  riGni  semi- -soan  and  tne  9tt6ct  of  its  imaoe  on  tne  isft 


semi -span  may  c-e  calculate 

the   s  Li  bscriot 


'eiationSa   i_pi 


sem 


sm 


SLiDscn  c 


represent   tne   righ" 


represer 


semi —span .   naki no   ti 


substitutions   for   the   symmetr 


e  i  e  m  e  n  t  s 


'en  nq  terms  =  equat  i  on s  !,  ** .  V4  )     ano 


m a v  c e  rewritten  a s 


.1/1) 


b&+bt$Lsintf  ij^fi 


I r, 


r,-.2 


=  ^^ .) 


4ir  Aft 


+ 


.(v-^if)       (vnp) 


rtedo 


EouatiGns  (7.2)  .   w.2-h-?  and  (7.25)  >  when   evaiaated 


pacn  control  do: 


t  r=  i=     =^  i=  m  ■ 


uuations       in       thp       sixnv — s-our       unknown 


rem  at  i  on  r  u.nc 


ne   solution   to   this 


quations  produces  une  desired  vaiu.es   of   the   circuiatxo- 

j.nct ion  at  the  mesh  points  over"  the  ri^ht  semi- span. 


BASIC  LIFT 


t  n  e  nasi 


t  orobl  8fn  a   pressure   distribution 


winq.   ieav_.no   the   correspond! nq   wi 


;1  ope  function  to  be  determined.  In  this  case,  equatio? 
6.22)  ,  rearranged,  gives  the  chordwise  derivatives  of  ths 
irculation  -function  at  each  mesh  noint  on  the  wi  nq  =  a; 
oil  o w s 


f> 


i  i  .  i  .-■  = 


AR 
8 


(-j)siaft  P(^>) 


Ht  une  leadzno  eone.  equation  i. 6 ■  i.b J  a.DDiie; 


AR 


r#,_<j»  =  f  c(^)PL(i) 


grp   pxnre s si o n  s 


i  nese 


Now ;  equations  £7. iO) f  (7. 14) ,  and  (7= 16) 
-For  \A    in  terms  of  the  unknown  circulation   values 
equations  are  repeated  below  with  known  values  collected  on 
the  left  side  of  each  equation. 


fy(j¥i>  +  (A0/4)  r^L(j>  =  (i/A^6nj,i} 


s: 


r0(j,i)    =    (l/A^)C<-l/2)P(j,i-l)    +    ii/2)V: 
ty<j,8>     =     (l/A0)C(-3/13>r<j,7)     -     (3/13)^ 


i  a  n  o 


26) 


bglu  res 


Now,     = 


CL!3C1  OH  5 


7)  .      ( 


=      J.   '_*  /      =  5SU   . 


r=!  i-nl 


i  c  L!  j.  a ts     ;_  n  c;    5D3 n  w  iss    q s?  i  v  a  1 1  v <=  s       c  ! 


<=  n-f  4-Ksa 


e  q  u  a  t  i  o  n  = 


'  j  i      1  — l  r~\  j— I  f       /  .-  — .     1  -t.  *#—  r— \.  *~\  t—\  i,i  > .-  *- •    j— .  t  .  •  »-i  I     r— i  ^— .  i-i-i      t— v  j— • 

.  —   .'  Ol  R  '_!  \    /     m    JL.  ■—!  .«  C>.  I       1=  !    I  O  *"»  r-.  1    5  *_'  W !    •    i  !    i    I  !=  »»1   HO 
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VIII.   SUMMAky 
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JNULUSi  UNs 


Hitnough  no  numerical  results  a.r&  yet  avail aoie,  there 
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calculation  method  presented  in  this  thesis.  A] I 
derivations  are  based  on  sound  mathemati ca3  procedures  and 
assumoti ons.   This  analysis,  alone-  with  the  parallel  effort 


nrpspntPd      i n 


ia.5      58'15tI90      Tine      1 


'   A     L_   J.    CT  . 


Durposi 


ie    c-^'^rs. 


research   et-fort,   which   was 


aoeouate  mathematical  model  ^-•D-r~    calculating  the  aerodvnam 
performance  of  a  wine  of  arbitrary   plan-form   in   invisci 


i n c omp  r e s s i d 1 e  +iow, 


J  .   L 


Lit;!   Uh  Kti~  tHhw. 


!_i  "="113.00  burf  acs  !  nsorv  "for-  winos  ow 


Arbitrary  i-'l  anforrn ,  Engineer's  Thesis,  Naval 
Post  or  adua  t  e  School  ,  Monterey  •  California,  1976. 


f-errin,  J.  J,  and  hjmirh,  L.  ;=>, 
Enoineers,  Prent i ce— Hal  1  ,  Inc=  ,  1979 


fter  od  vnami  c  s  -for 


M.  and  Chow.  C 


r  Dunaati  ons  o- 


4, 


fierDdynaiDics:  Bases  of  nerooynamic  Design, 
Third  Edition,  Wiley,  1976= 

Barber,  E.  M. ,  A  Method  to  Calculate  the  Paramsti 
of  Wings  of  Arbitrary  Planforro,  M.S.  Thesis,  Na.v, 
Postgraduate  School,  Monterey.  California.  1934. 


58 


INITIAL  DlbTRiHUTIUN  Libj 


Uetenss  sechnicai  Information  Center 

Cameron  Station 

Alexandria,  Virginia   22314 


No 


Library*  Code  U142 
Naval  Postgraduate  School 
onterey.  California   93943 


M 


Professor  T.H,  Gawain,  Code  67  bn 
Department  of  Aeronautics 
Naval  Post Graduate  School 
Monterey?  California   93943 

Chairman,  Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California   93943 

LT  E.M,  Barber 

FAIRECONRON  THREE  (VQ-3) 

NAS  Barbers  Point,  HI   96862 


LCDR  Christopher  S,  Jones 
jsj  A  V  A I R  S  V  S  C  0  M  ^  A I R — 5  3  6  > 


wasni no con  = 


n  r 


20361 


59 


Thesis 

J68U5 

c.l 


Jones 

A  computational  me- 
thod for  wings  of  arbi- 
trary planform. 


Thesis 

J68U5   Jones 

c.l         A  computational  me- 
thod for  wings  of  arbi- 
trary planform. 


